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The ^H(p, n) charge exchange reaction at Tp = 790 MeV is used to study the A(1232)-nucleon 
(AN) interaction in the A resonance excitation energy region. For the AA'^ potential, a meson 
exchange model is adopted where tt, p, ui, and a meson exchanges are taken into account. The 
deuteron disintegration below and above pion threshold is calculated using a coupled channel ap- 
proach. Various observables, such as the inclusive cross section, the quasifree A decay, the coherent 
pion production, and the two-nucleon breakup are considered. It is shown that these observables are 
influenced by the dynamical treatment of the A degrees of freedom. Of special interest is the coher- 
' ent pion decay of the A resonance which is studied by means of the exclusive reaction ^H(p, n7r"'")^H. 

I Both the peak energy and the magnitude of the coherent pion production cross section depend very 

■ sensitively on the strength of the AA'^ potential. The coherent pious have a peak energy of uj — 300 

MeV and a strongly forward peaked angular distribution. 
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I. INTRODUCTION 



<N . . . , 

^ ' In recent years, inclusive and exclusive {p,n) and { He,t) charge exchange reactions at intermediate energies have 

On , been proven to be excellent probes for investigating the A dynamics in nuclei. The most important observation 
of the inclusive charge exchange reactions is the downward energy shift of the A resonance peak position by ~ 70 
MeV in nuclear targets (with mass number A > 10) as compared to the proton target 13-0]. Microscopic A-hole 
model calculations show that a lar ge p art of this shift is caused by the attractive A-nucleon [AN) interaction in the 
spin-longitudinal (LO) channel |]6| |lO|]. A direct signature of this interaction is provided by the measurement of the 
^\ I coherent pion decay spectrum where the pion of the A resonance decay is measured in coincidence with the ejectile 
while the target nucleus is left in its ground state |pd]-|l6|. These coherent pions couple strongly to the LO channel. 
Due to the attractive AA^ potential in this channel, their energy spectrum is substantially shifted downwards relative 
to the A resonance peak position of the inclusive reaction. 
O , In this paper we study the A excitation in the deuteron target using the ^H(p, n) charge exchange reaction at 
2 rp=790 MeV. The deuteron has the advantage that its wave function is well known and that the Fermi motion of the 
AA^ system can be treated properly. Therefore the effects of the AN interaction can be studied in a more direct way 
than in heavier nuclei. In the past most of the studies on the deuteron have been carried out with electromagnetic 
and hadronic probes, such as in photon-deuteron {'jd) and pion-deuteron (jrd) scattering, leading to NN, nd. and 
ttNN final reaction channels. In these reactions the intermediate AA^ interaction plays an important role [0-g3|. 
Ci I The photon excites the A dominantly with spin-transverse (TR) coupling, i. e. by the transition operator S x qT 
{S and T are the spin and isospin transition operators, respectively), while the pion excites it with spin-longitudinal 
(S ■ qT) coupling. Both couplings are orthogonal to each other and therefore give different information on the AA'^ 
interaction. In the charge exchange reactions the target is exposed to the virtual tt and p meson fields produced by 
the {p,n) projectile-ejectile system. The pion-like interaction excites the LO response function {S ■ qT coupling) 
of the target while the p-meson-like interaction excites the TR response function {S x qT coupling). Due to the 
kinematics, the virtual meson fields obey the energy-momentum relation uj < q and thus explore the LO and TR 
response functions in an (w,(f) region which is inaccessible to real pion and real photon scattering. 

In the interpretation of the (7^) and (ttc?) scattering reactions, various theoretical models have been applied to 
treat the AA^ dynamics. Among these models are the coupled channel approach |24|~[29[| and the three-body Faddeev 
treatment of the ttNN system ||3^,^ . In the present paper we make use of a coupled channel approach to describe 
the AA^ system in a non-relativistic framework. We set up a system of coupled equations for the AA^ wave function in 
configuration space. We then apply the Lanczos method for solving the equations. The AiV potential is constructed 
within a meson-exchange model. The exchanged mesons taken into account are the pion (tt), the rho (p), the omega 
(w), and the sigma (a). 

The aim of the present paper is to show that the {p, n) reaction at forward scattering angles is an interesting tool 
to study the AA^ interaction in the deuteron. The advantage of the (p, n) reaction over other probes is twofold: 
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First, at forward scattering angles the quasi-elastic peak cross section and the A resonance peak cross section are 
energetically well separated. The A resonance cross section is very large as compared to the quasi-elastic cross section 
providing a well-defined A resonance peak in the spectrum. Second, the A resonance cross section involves both a 
LO and TR component (with a ratio of LO/TR =1/2) and thus allows to examine the complete spin structure of 
the AA^ interaction. The LO excitation is of special interest for the coherent pion decay of the AA^ system. Due to 
the LO spin structure in both the excitation and de-excitation process the cross section for coherent pion production 
becomes relative large. Because the AiV interaction is most attractive in the LO channel, the coherent pion spectrum 
is expected to clearly exhibit the effects of the AiV potential. It turns out, indeed, that both the peak energy and the 
magnitude of the coherent pion production cross section depend quite sensitively on the strength of this potential. 
Other observables, such as the inclusive cross section, the quasifree A decay cross section, and the two-nucleon 
breakup, are also shown to be influenced by the dynamical treatment of the A degrees of freedom. 

The organization of the present paper is at follows. In Sec. 2 we give a detailed account of the formulation and 
methods of calculation used in the analysis of the data. First we present the coupled channel approach and show how 
to calculate the correlated AA^ wave function in a very efficient way with the Lanczos method. Then we split the 
inclusive cross section into its various partial cross sections, as there are the contributions of coherent pion production, 
quasifree A decay, and two-nucleon breakup. In Sec. 3 we discuss the parameters used in our model. We present the 
results of the cross section calculations and compare them to experimental data p2| , ^ . Finally, in Sec. 4 we give a 
summary and conclusions. 



II. THEORY 



We are interested here in the calculation of the inclusive and exclusive cross sections for the ^H(p, n) charge exchange 
reaction in the A resonance energy region. Since we shall deal with high projectile energies {Tp = 790 MeV), we assume 
that the cross sections can be calculated within the distorted wave impulse approximation (DWIA). Because of the 
transparency of the deuteron target it is sufflcient to calculate the distortion effects within the eikonal approximation. 



A. Inclusive cross section 



We start our formulation by writing down the formula for the inclusive ^H(p, ti) cross section. Using relativistic 
kinematics, this cross section is given as 

2EpEa MpMn d^pn , ^ p p ^ i |2 n\ 

d<^= I , T^7^ r9^^3 9F 2^{^TT)5(Ep + Ed-En-Ef)\Tf, \ . (1) 

l\{s,M^,Mj) ^" (2n)^2En^ 



Here the indices p{n) and d refer to the proton projectile (neutron ejectile) and the deuteron target, respectively. Ef 
is the total energy of the final state on the target side, i.e. of all outgoing particles except the neutron ejectile. The 
cm. momentum of the final state is fixed by momentum and energy conservation while an integration over the relative 
momenta has to be performed. In addition, an average over the initial spin orientations and a sum over the final spin 
orientations of both the projectile and the target spin states are taken. The transition amplitude will be evaluated 
in the Breit frame (BF) of the target system, where the deuteron state is well described by the usual nonrelativistic 
wave function. 

Since we want to compare our results with experimental data where the energy and the scattering angle of the 
outgoing neutron have been measured, we have to calculate the double differential cross section d^a/dE^diln in the 
laboratory frame (lab) . From Eq. (0) we obtain 

d^a X'""^ p\f / d^a 



dEndQn J Pn XdEndQ 

In the following, cross section results will always be given in the laboratory system while every other unlabeled 
quantity refers to the Breit Frame. 
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B. The uncorrelated source function 



In this section we discuss the excitation processes which contribute to ^H(p, n) in impulse approximation. In our 
model, there are two relevant Feynman diagrams for this case which are shown in Fig. ^ They represent the nucleon 
excitation (a) and the A excitation (b) from the deuteron target. In first order DWIA the transition amplitude Mfi 
corresponding to the graphs in Fig. m is given as 

Mfi = X! / ^^^0 4>i^^*iPniro) (Xn I (V'afc I toj{uj,fQ - Tj) \ V^d) I Xp) ^p"^^ (Pp, ^i)- (3) 

Here rj denotes the coordinates of the projectile (j = 0) and the target nucleons {j — 1,2), respectively, which are 

measured relative to the center of mass of the target. cj)p and cjrn are the distorted wave functions of projectile 
and ejectile in the initial and final channels. The spin-isospin part of the projectile (ejectile) wave function is denoted 
as I Xp) {{xn I )• The initial and final states of the target are | ijjd) and {ipab \^ where ah refers to either the NN or 
the AA^ system. Note that Mfi describes only the excitation process but not the de-excitation, e.g. the free decay of 
the A. Therefore M fi has to be distinguished from the transition amplitude Tfi for the complete reaction. 

The toj in Eq. (||) is the effective interaction between the projectile nucleon and the target nucleon j. A sum 
over the two target nucleons has to be performed. The toj is represented by the free nucleon-nucleon tpf^^^N matrix 
in the case of the nucleon excitation (Fig. ^) while it is approximated by the free A^A^ NA transition operator 
tNN,NA for the A excitation (Fig. |l|b). The specific form of these interactions is given in App. They provide 
a spin-longitudinal (pion-like) and a spin-transverse (p-meson-like) excitation component. The excitation strength 
parameters are fitted to reproduce the experimental cross section and the spin observables for (p, n) reactions on a 
nucleon target. 

Since the effective toj is adjusted to experiment, we implicitly account for knockout exchange effects between the 
projectile and the target nucleons. Note, however, that we neglect the projectile A excitation process where the 
detected neutron comes from the decay of the A. For forward neutron angles, the cross section contribution of 
projectile excitation gives only a small correction to the dominant target excitation process. In Refs. [ p4[|35| ], the 
projectile excitation is found to be suppressed by a factor of ~ 10 in the A energy region. 

We now take advantage of the fact that at high incident energies and large momentum transfers both interactions 
turn out to be rather short ranged, i.e. very weakly dependent on the four-momentum transfer (w, q) = {Ep — E^Pp — 
Pn). Therefore they can be well approximated by local operators in r-space of essentially 5-function form, that is 

tf)]{uj, ro - r-j) = -f^;^ j '^^l' exp[ig" • (rp - fj)] toj{cJ, q') w taj{uj, q) (5^(ro - r^). (4) 

By inserting Eq. (||) into Eq. (^) we find for the transition amplitude 

Mf, = {i^ab I /3 I Vd), (5) 
where the hadronic transition operator p is defined as 

P= ^{Xn \ toj{uj,q) XDw{Pp,Pn,rj) I Xp), (6) 

with 

Xdw {Pp ,Pn,rj) ^ {pn , rj )<l>p{Pp,rj). (7) 

The product wave function Xjjw includes the projectile and ejectile distortion effects in the reaction. In the eikonal 
approximation, Eq. (|^) is reduced to XDwi<l,^j) — ^dw{q) exp {iq - fj ), where Njywiq) can be calculated from the 
distorting potential. Thus we get from Eq. ([g) 

P = Nowiq) iXn I toi{uj,q)exp{iq- r/2) + to2{uj,q)exp{-iq- f/2) \ Xp), (8) 

where r = ri — ^2 is the relative coordinate between the two target nucleons. 

Following Ref. we can now introduce source functions for the excitation of either a NN or a AN system, 
respectively, by defining 

I PN{f)) = Nowiq) {Xn I C'^'^^iuJ, q) exp (zg • f/2) \ Xp) \ Mr)) + (1 2), (9a) 
I PAif)) = Nowiq) iXn I C^'^^ito, q) exp {iq ■ f/2) \ Xp) \ Mr)) + {l^2). (9b) 
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The source functions represent the doorway states excited initially by the external (p, n) charge exchange field. We 
call them uncorrelated since they do not include the "final state" interactions within the excited NN or A7V system. 
If the final state interactions are neglected, the second target nucleon is just a spectator which does not take part in 
the interaction with the projectile ("spectator approximation"). The full dynamical treatment of the NN and A7V 
systems requires the calculation of correlated continuum wave functions which we discuss in the next section. 



C. Inclusion of the AN interaction and calculation of the correlated wave function 

The inclusion of final state interactions in the description of the ^H(p, n) reaction leads to a coupled channel 
problem. In the energy regime considered here, up to four particles appear in the final state because of the pion 
production. Therefore the full coupled channel problem is complicated and can only be solved approximately. In our 
model, we will treat the A as a quasi-particle with a given intrinsic decay width and mass. The configuration space Ti. 
is build up from NN, AN, and irNN sectors, Ti = TCnn ®'Han (BH-nNN- The corresponding projection operators are 
denoted as Pn, Pa, and Pq, respectively, and we will use for any operator Q the obvious notation ^Ian = Pa^Pn, 
^QA = Pq^Pa, etc. 

In this formulation the source functions of Eqs. ( ^J9b| ) appear to be projections of a unique source function | p), 
i.e. 

Pn\p)^\pn), Pa|p)=|pa). (10) 

Since we restrict ourselves to DWIA there is no source function for the ttNN sector, i.e. Pq \ p) ~ Q. We now define 
the correlated continuum wave function | ip) as 

I ^) = z, l^- \P) = G\P). (11) 
h — H + le 

The Hamiltonian H = H'^ + V contains the free energy of the system as well as the interaction V. The full 
propagator G is connected to the unperturbed Greensfunction = {E — H°)~^ by the equation 

G = G° + G°VG. (12) 

The interaction V introduces the correlation effects into the wave function | We will demonstrate next how to 
solve Eq. (|l^) in a very efficient way by using a few simplifying but physically reasonable approximations. 

First we make use of the fact that pion production in intermediate energy charge exchange reactions is dominated 
by the A resonance excitation. We assume that Vnq = Vqn = 0, i.e. there is no direct coupling between the NN 
and the ttNN sector. This means that in our model pion production proceeds exclusively via an intermediate AA^ 
system while the nucleon pole terms are neglected. In the A excitation energy regime this approximation should be 
very good. 

The correlated AiV wave function | iPa) follows from Eq. ( pi] ) by projecting on the AN channel, i.e. 
I iPa) — Pa I V') • the A energy region, | ^Pa) can be well approximated by 

\M~Gaa\pa). (13) 

Here we neglected the contribution Gan \ Pn) from the NN source function which is substantially suppressed in the 
A energy region. The propagator Gaa of the AA^ system is determined from Eq. (|lj) which yields 

Gaa = G^A + G^a^aaGaa + G^^VaqGqa + G^^VanGna- (14) 

In this equation, the first and second term on the r.h.s. represent the free propagator and the interaction term 
in the AN channel, respectively. The third term G^^VaqGqa refiects the emission and re-absorption of pions 
(A — > ttN — > A) which leads to the energy dependent decay width rA(sA)- However, since we treat the A as a 
quasi-particle with an intrinsic width and the physical mass Ma = 1232 MeV, this "self-dressing" process is already 
effectively included in the free propagator . Therefore the third term of Eq. ( |l4|) has to be dropped. In addition, 
we also drop the fourth term G^^VanGna which generates NN box contributions to the AN potential. We do not 
expect those contributions to influence the AiV interaction significantly since it seems to be much more likely that 
the A is produced directly in the interaction with the projectile than later on in the final state interaction of two 
outgoing nuclcons. Hence we are left with 

Gaa ~ G'^aa + G^^VaaGaa (15) 
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for the AA'^ propagator. This approximation has the major advantage that the AA'^ channel effectively decouples from 
the NN channel and therefore can be solved separately. 

If we project now Eq. (Ilj) to the possible final channels NN or t:NN and make use of Eq. (Esf), we find 



with 



Ppf \ip) = Gnn I Pn) + GnnVnaGaa \ Pa), (16a) 
Pq\^) = GqqVqaGaa \ pa), (16b) 



Gnn — G%j^ + G%j^VnnGnn, (17a) 
Gqq = ® Gnn- (17b) 

Transitions from the AA^ to the NN or ttNN channel are accounted for by the corresponding transition potentials 
Vna and Vqa in Eq. (|l^). They appear only in first order. The potential Vnn in Eq. (|l7a|) describes the NN final 
state interaction which includes also AA^ box contributions. We adopt the Paris potential i36| , |3^ for Vnn- All effects 
of the A^A^ po tenti al can be included in the A^A^ wavefunction {ipNN \ for the outgoing nucleons. Furthermore, we 



assume in Eq. ( 17b ) that the outgoing pion is not distorted from the remaining A^A^ system so that it can be described 
by a plane wave. 

In the remainder of this section, we focus on the AA^ subsystem and the method for solving | -^a) ~ Gaa | Pa) ■ 
After having separated the cm. motion, the full propagator for the relative AA^ wavefunction is given by 



with the excitation energy 



and the relative kinetic energy 



Gaa = ■ ; — , (18) 

CA + 2rA(sA) - Ta - Vaa 



In the spirit of the model, the propagator contains the energy dependent width rA(sA) of the A resonance. The 
invariant mass sa is fixed by conservation of four-momentum at the production vertex. For given projectile kinematics 
we obtain 

SA = (M + C^lab)'-g;ab- (21) 

Here to is the energy transfer and q the momentum transfer to the deuteron target. Note that the invariant A mass 
is calculated by assuming the target nucleon to be at rest in the laboratory frame (frozen approximation) . 
In order to solve | ipA) ~ Gaa \ Pa), wc transform this equation into an equivalent integral equation |10| 

|Aa)=|pa)+1^aaG1a|Aa), (22) 

so that 

|^a) = G^a I Aa). (23) 

Eq. ( ^^ is now reduced to a set of coupled channel equations for radial wave functions in the following way. First 
we expand | pa) and | Aa) in terms of partial waves, i.e. 



Pa)= l:PSL.jMAr)\iSL)JMj)\n), (24) 

SLJMj 

Aa)= -^SLJMAr)\{SL)JMj)\n), (25) 



SLJMj ^ 
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where S, L, J, Mj denote the spin, orbital and total angular momentum quantum numbers of the AA'^ system. Due 
to the special isospin structure of the excitation process, the isospin channel is always | TMt) =| 11) • After insertion 
of Eqs. (HH) into Eq. (|2|) and projection on {{S'L')J'M'j \ we obtain 

^A„(r) = ip„(r)+^ I dr'r'2K„'(r)Gi^(r,r')^A„,(r'), (26) 

n' 

where we used the index n as a shorthand notation for {SLJAIj}. Eq. (E9) may easily be written as a matrix 



equation. The radial source functions p„(r) can be calculated in a straightforward manner and are given in App. B 1 
The matrix elements Vnn' (f) of the AA^ potential, which is discussed in the next section, arc given explicitly in App, 
B2. Note that the operation of onto A„'(r') involves a radial integration besides the matrix multiplication. 



The merit of solving for | Aa) first lies in the fact that the corresponding radial functions A„(r) are localized. This 
fact makes it possible to apply the Lanczos method which is described in Ref. It turns out that the Lanczos 

method allows us to solve Eq. ( p6| ) in a very efficient way. Once | Aa) is known, it is easy to calculate | ipA) from Eq. 
( p3| ) as well. 

III. THE MESON EXCHANGE MODEL FOR THE AiV INTERACTION 

Similar to the NN interaction, the AA^ interaction can be constructed within a meson exchange model. In the 
present work we use tt, p, lu, and a exchange. While the to and a mesons contribute only to the direct term shown in 
Fig. ^(a), the tt and p mesons may also induce spin-isospin-flip transitions which lead to the exchange term of Fig. 
|2|(b ). We will discuss the latter term first and start from the interaction Lagrangians C-^na and CpNA given in Ref. 
|]38[ . In the non-relativistic reduction we obtain 

v:Ak) = % T\ . ,f\:f^v^^. + (1 - 2), (27) 



V£.ik) = % • + (1 ^ 2). (28) 



^0 - n, -m-+te 

The corresponding direct interaction terms follow from these potentials by replacing the transition matrices S, 
{TjT'^) with the spin (isospin) 3/2 matrix S (9) and the spin (isospin) 1/2 matrix a (r), respectively. This yields 

V,lik) = ^^^^ e, . (^^^)(^ -^) + (1 ^ 2), (29) 

Vi,ik) = If^^ e. . r1 ^I'^'^^Y'^ + (1 - 2). (30) 
'Tip k^ — k'^ — rrip + le 

For explicit calculation of these potentials we need to know the appropriate coupling constants. The ttNN and 
pNN couplings are well known experimentally from an analysis of NN scattering data pq] . They may be related to 
the other couplings by 

f-rr/pNA — '^f-n/pNNi (31a) 
/tt/pAA = -^Iir/pNN- (31b) 



Eq. (|31m fo llows from the static quark model |39[|. For Eq. (31a) we decided to use the widely accepted Chew-Low 
relation |40|| instead because the quark model prediction for /t^na is too small as compared to the experimental value 
from the A decay. 

Additional contributions to the direct part of the AiV interaction emerge from the uj and a exchange. We assume the 
couplings ljAA and crAA to be the same as for ojNN and aNN, respectively. As before, we choose the non-relativistic 
limit of the interaction which follows from C^^nn and C^nn [pq]. The resulting potential is 



6 



(32) 



The full potential for the AA^ interaction is the sum of the different contributions from Eqs. (27 - ^) and (|3^), 
hence 

VAAik) = Kl^^ik) + V^+'ik) + V-+^{k). (33) 

Note that the tt (p) contribution covers the LO (TR) spin-isospin part of the potential while the lo -\- a contribution 
is spin-isospin independent. Furthermore, we point out that we use monopole form factors of the type 

F{e) ^ { (34) 

at all vertices. In the formulas given here, the form factors are always suppressed in order to simplify the notation. 

Since the AiV potential is required in r-space, a Fourier transformation of VAA{k) has to be performed. This 
involves an integration over all possible three momentum transfers k. On the other hand, the energy transfer fco is a 
free parameter which has to be fixed within the kinematical allowed region. Because the A resonance is assumed to 
keep its invariant mass during the propagation, wc choose fco = in the direct interaction terms and fco = w in the 
exchange terms. This choice makes the potential energy-dependent, as it should be. 

The AiV NN transition potential Vat a is constructed in complete analogy from Eqs. (^, |2^) by replacing 
(TT) with (T (r). Of course, only the tt and the p meson contribute to Vata- For the energy transfer in this case, we 
adopt the choice /co — ^l'^. of Ref. [0. As far as the spin-longitudinal part of Vna is concerned, we use an additional 
zero-range interaction 

V^^{k) = gj^J-^Il^hHK fi . r1 {S^ ■ k)iaV ■ k), (35) 
mi 

with the Landau-Migdal parameter qna = 1/3 [p^ . This additional contribution cancels the non-physical zero-range 
part of the tt exchange potential. In principle one should introduce a corresponding Landau-Migdal term for Vaa with 
a strength parameter gAA e^s well. However, due to the presence of not only exchange but also direct contributions to 
the potential, the zero-range parts cancel anyway and (?aa turns out to have no influence on the results. Therefore 
we decided to do without and chose ^aa = 0. 

In Tab. || an overview of all the meson parameters used for the A7V potentials is given. We made sure that 
with these parameters, our model consistently reproduces experimental results for pion absorption on the deuteron. 



IV. DECOMPOSITION OF THE INCLUSIVE CROSS SECTION 



We will now decompose the inclusive ^H(p, n) cross section into partial cross sections corresponding to different 
physical processes. These processes are schematically represented by the diagrams of Fig. |(a)-(d). For each process, 
only the lowest-order diagram is shown. We distinguish between quasi-elastic scattering (a), p-wave rescattering (b), 
coherent pion production (c), and quasifree A decay (d). 

Quasi-elastic scattering and p-wave rescattering both result in a two-proton {2p) final state on the target side. 
Their transition amplitudes interfere coherently with each other. While the quasi-elastic scattering directly leads 
to the 2p system, the p-wave rescattering contribution involves an intermediate AA'^ state and arises due to the 
AN — > NN transition potential Vna, see Fig. H(a,b). Since in both cases the deuteron is simply broken up into two 
protons, we will refer to this reaction channel as to the breakup channel (BU). 

At higher energy transfers, the nNN final channel is the most important one. We remind the reader that we assume 
the pion production to proceed always through A resonance excitation. If the two nucleons in the final state form a 
bound deuteron state, we speak of coherent pion production (CP), see Fig. |3|(c). In this case, the outgoing pion is 
naturally a tt'^. On the other hand, if the two outgoing nucleons are unbound, we speak of quasifree A decay (QF), 
see Fig. |3|(d). There are two possible final isospin configurations for this process, namely ■n~^np and ir^pp. 

According to these different reaction mechanisms, the inclusive cross section of Eq. ^ may now be split into its 
various components. We express the squared transition matrix as the following sum, 

J2 ^{^ + Ed- Ef) I Tf, |2 = i {Sbv + Sep + 5qf) . (36) 
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Here the strength functions S impUcitly include the integration over free relative momenta as well as the summation 
over all spin confi gurations, and the factor g arises from the average over the initial spin configurations. 
Using Eq. ( |l6aD , the strength function for the deuteron breakup channel BU is found to be 

/,3 
— ^ 5{uj + Ed- E2p) I {^NN I pn) + {ipNN I Vna I V-a) I', (37) 

The first matrix element in the sum of the r.h.s. describes the quasi-elastic scattering, and the second matrix element 
describes the p-wave rescattering. p denotes the relative momentum of the two outgoing protons. Their final state 
interaction Vnn is included in {ipNN \ while the correlated source function | '0a) accounts for the AiV interaction 
Vaa- 

In order to calculate the matrix element for coherent pion production, we define the operator 

= e-'«-' -^^/^ Ti, . + (1 ^ 2) , (38) 

rrirr 

where is the relative pion-nucleon momentum in the A rest frame. The operator describes the decay of the A 
into a real pion and a nucleon and contains also the plane wave of the outgoing pion. Applying /tt to the AA^ wave 
function and projecting onto the final deuteron state yields 



ScP^J2 j J^^^Ji^ + Eci-E^-E',) \{^a\h\^^A)? (39) 



I Tn 

S..Sf' 

for the CP strength function. Eq. ( ^ ) may also be rewritten as 

Sep ^y2{^A\ -MVaa) I 0a), (40) 

where we used the fact that Ed ~ E'^ holds in the Breit frame. Hence the coherent pion production amplitude is 
directly connected to the imaginary part of the AN interaction. We remark that this is only true if Im(VAA) comes 
exclusively from the tt exchange interaction, as it is the case in the A energy region. 

From Eq. ( |39| ) we may easily calculate the angular distribution of the coherent pions, namely the triple differential 
cross section dracp/dEnd^lndflTr, by omitting the integration over the pion momentum. The result is 

d^acP M,M„p'-^Ed ql dq, f , 

m I u I V'a) I (41) 



dEndnndn^J (27r)5 p^^^ Md 2E^ dEf 

with the same notation for the phase space factor as in Eq. (||). 

In analogy to Eq. (^0|), the pion production from quasifree A decay is ascribed to the imaginary part of the A 
"self-dressing" , i.e. to the intrinsic A decay width. As a consequence, the QF strength function is given as 

^QF = E(V^a|-^^|0a). (42) 

The differential QF cross section could also be calculated from the matrix element | {ipNN \ fv \ V'a) Eq- (^) would 
then be recovered after integration over the pion momentum, as in Eq. (^9|), and over the relative NN momentum as 
well. 

To summarize we have decomposed the inclusive ^H(p, n) cross section into three parts: the quasifree A decay 
(QF), the coherent pion production (CP) and the two-nucleon breakup of the deuteron (BU = quasi-elastic plus 
p-wave rescattering) . Of course there are other possible contributions to the inclusive cross section, but in the energy 
region under consideration, they should be less important than the ones mentioned and are therefore neglected in the 
present paper. 



V. RESULTS AND DISCUSSION 
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A. Input parameters 



With the formahsm described in the previous sections, we have calculated energy spectra at various scattering 
angles On for the ^H(p, n) reaction at Tp — 790 MeV. The deuteron wave function and the two-nucleon wave function 
were generated from the Paris potential [^,|^. A constant distorsion factor of N'^^ = 0.9 was applied p2| , ^ . 
For the AA^ potential we used coupling constants and cutoff parameters as given in Table |. The free A mass is 
Ma = 1232 MeV, and the free decay width rA(sA) was parameterized in the usual form |^^. 

The strength parameters in the effective parameterizations of ^ATAr.TVAf and Inn^na (cp. Sec. [IB, App. ^) have 



been fitted to reproduce experimental results for the nucleon target. In order to demonstrate the quality of this fit, 
we present in Fig. ^ results for the -'^Hfn, n)A++ cross section at Tp = 790 MeV. The theoretical calculations are 
compared with the experimental data ||32|] at various scattering angles On of the outgoing neutron. Both shape and 
magnitude are reproduced very well. We remark that also the LO/TR ratio of 1/2 used in the parameterization of 
tNN.NA has been observed experimentally [H,E5|. 



B. Inclusive spectra of the •^H(p, n) reaction 

In Fig. H we show the calculated inclusive cross section for the On = 0° spectrum of the ^H(p, n) reaction in 
comparison with the experimental data [3^ . The different contributions to the inclusive cross section are also shown 
separately. The possible reaction channels are the quasifree A decay (QF), the coherent pion production (CP), and 
the deuteron breakup (BU) due to quasi-elastic scattering and p-wave rescattering. 

At energies wiab < 50, the spectrum exhibits a prominent peak which arises dominantly from quasi-elastic scattering. 
The small width of the peak reflects the Fermi motion of the nucleons in the deuteron target. The peak is described 
well by our calculations. 

In the energy region of the A resonance peak position at wiab — 340 MeV, the theoretical result is also in good 
agreement with the data. The cross section in this region is mainly due to the quasifree decay of the A. The cross 
section contribution due to coherent pion production is by a factor of « 5 smaller as compared to the QF contribution. 

At high energy transfers wiab > 450 MeV, our calculation underestimates the experimental cross section. The miss- 
ing cross section can be ascribed to the excitation of higher nucleon resonances, e.g. N*{1A40). Those configurations 
are not included in our model space, and therefore the underestimate of the cross section is not surprising. 

The theoretical calculations also underestimate the data in the so called dip-region between the quasi-elastic peak 
and the A peak, i.e. in the energy region 50 MeV < wiab < 250 MeV. The experimental cross section here is believed to 
result mainly from two-body exchange currents in the target. Only the baryonic exchange currents connected with the 
A are accounted for in our calculations. They are included in the AA^ interaction Vaa and in the transition potential 
V/vA which gives rise to the p-wave rescattering contribution. In Fig. ^, our full model calculation is compared with a 
calculation obtained within the "spectator approximation". The latter one corresponds to the case Vaa = Vna — 0. 
It can be seen that the inclusion of the exchange currents significantly improves the theoretical description of the 
experimental data. Nevertheless, we still underestimate the cross section in the dip-region by a factor of about 1.5. 
This may be due to the neglect of purely mesonic exchange currents, e.g. the s-wave rescattering of pions. At the 
low energy side of the A peak, also relativistic corrections in the dynamical treatment of the AA^ system could be of 
some importance. 

In the following we shall demonstrate that almost all of the cross section in the dip-region is spin-transverse (TR) 
while the spin-longitudinal (LO) cross section is relatively small. This statement can be tested by measuring spin 
observables that allow for a separation of the inclusive cross section into its LO and TR components This 
separation is shown in Fig. ^ together with the theoretical results. The measured LO cross section is reproduced very 
well by our calculation. On the other hand, the measured TR cross section is both shifted and enhanced with respect 
to the calculation. This apparent shift clearly indicates that our calculations produce not enough TR cross section in 
the dip-region. A similar effect is observed in inelastic electron scattering off nuclei |^,^. In the (e, e') reactions, 
the virtual photon also probes the TR response of the target. In scattering off nuclei, the experimental (e, e') cross 
section in the dip-region is considerably larger than expected. The additional cross section cannot be explained by 
one-body currents and has to be ascribed to two-body currents Q . Such rescattering effects have also proven to be 
quite important for coherent pion photoproduction on the deutero n pl| ]. In hadronic interactions, the p meson should 
take over the role of the photon. In a recent paper of C.-Y. Lee |)49|, the author found indeed a large contribution 
to the ^II(p, n) inclusive cross section from the exchange of a p meson which couples to a "pion-in-flight" in the 
deuteron. We could not reproduce this result. Since there are many other possible meson exchange currents which 
have not been considered for the ^H(p, n) reaction so far, the dip-region puzzle remains a problem to be solved. 



9 



The theoretical description of the experimental ^H(p, n) spectra at higher scattering angles is of the same quality 
as for zero degree scattering. This can be seen from Fig. |^ where we compare the theoretical ^H(p, n) cross section 
with data at On = 7.5° and 0„ — 15°, respectively. At the low energy side of the A, the theoretical curves show again 
a characteristic underestimate of the data. Apart of this effect both the shape and the magnitude of the spectra are 
described well. Our calculations correctly reproduce the change of the quasi-elastic and the A peak position with the 
scattering angle. 



C. Influence of the AA'^ interaction on the exclusive spectra 



We now discuss the influence of the AA^ interaction on the three partial contributions (QF, CP, BU) to the cross 
section. As we will demonstrate, the coherent pion production is most sensitive to Vaa- This sensitivity is a combined 
effect of the attractive pion exchange, the spin-longitudinal ttA'^A coupling and the structure of the deuteron wave 
function. In the following discussion, we will focus on the exclusive ^H(p, ?i7r+)^H spectra first. 

In Fig. H, theoretical results for the exclusive ^II{p, n7r"'")^II cross section are shown. The calculation without 
inclusion of the AA^ interaction (i.e. T^aa = 0) is compared with a calculation where only the tt meson contribution 
to Vaa was taken into account. The effects of the direct and the exchange part of the tt mediated interaction are also 
examined separately. Obviously both parts are of equal importance which may be surprising at first regard. It is true 
that because of the smallness of the ttAA coupling constant, the direct part should be suppressed by a factor of « 20 
as compared to the exchange part, but in the T — \ channel of the AA^ system which is relevant here, the spin-isospin 
matrix elements turn out to be larger by approximately the same amount. Therefore it is essential to include both 
the direct and the exchange contribution in the AA^ interaction. In the case of the tt exchange, the interference effect 
of both contributions leads to a very strong attraction between the A and the second target nucleon. The result is 
a shift of the cross section peak position downwards in energy (by ~ 60 MeV, cp. Fig. ||). Most of the attraction is 
caused by the tensor part of the tt whereas the central part is less important for the peak position but influences the 
overall magnitude of the result. 



By including not only the tt, but also the p, w, and a mesons in VaAj we obtain the results presented in Fig. 10. It 
can be seen that both the peak position and the magnitude of the exclusive ^H(p, n7r+)^H cross section depend quite 
sensitively on the specific form of the AA potential. In comparison to the result with only tt exchange, the inclusion 
of TT and p in the AA interaction leads to a less attractive potential and hence to a smaller shift of the peak position. 
The reason for this behavior is the partial cancelation of the tt and p tensor forces which have opposite sign. For 
the direct part of the interaction, the lo and the a meson cause an additional short range repulsion and a medium 
range attraction, respectively, which leads to an additional enhancement of the cross section. The final result for the 
^II(p, n7r+)^II spectrum (i.e. the result with inclusion oin + p + uj + am Vaa) is quite different from the spectator 
approximation (i.e. the result with Vaa = 0). Most remarkable is the shift of the peak position downwards in energy 
by about 30 MeV. As discussed above, the main reason for this shift is the attractive pion exchange, i.e. the attractive 
LO part of the AA^ potential. The amount of the peak shift is directly related to the interaction strength of Vaa- 

This result attracts further interest because there is no comparable peak shift due to Vaa for the quasifree A decay 
cross section and also not for the deuteron breakup. Fig. |ll] shows the results with and without inclusion of Vaa 
for the coherent pion production (CP), the quasifree A decay (QF) and the deuteron breakup (BU). In both the QF 
and the BU channel, the cross section is slightly enlarged and broadened due to the AA^ interaction, but there is no 
significant shift of the peak position. 

In order to understand this different behavior, we perform a multipole decomposition of the partial cross sections 
and split them into two contributions corresponding to unnatural parity (UP) states (J^ = 0~, 1"*", 2~, . . .) of the 
AN system and to natural parity (NP) states {J^ = 1~, 2+, . . .), respectively. For the three reaction channels under 
consideration, this decomposition is shown in Fig. |l^. For each case, the full calculation is compared with the spectator 
approximation. One recognizes that the cross section contributions of the UP states are always lowered in excitation 
energy by the AA^ interaction while the NP states are not. This is explained by the fact that there is a strong coupling 
of the pion to the UP (pion-like) states but just a weak coupling to the NP states. Therefore only the UP states are 
substantially influenced by the attraction of the LO part of Vaa ■ For the coherent pion production, the spectrum is 
clearly dominated by the UP states. This is an effect of the LO spin structure of the de-excitation process and of the 
deuteron wave function (which selects spin S — 1 and orbital momentum L = 0, 2) in the final state. 

For the quasifree A decay, the final deuteron state is replaced by a final unbound A^A^ state. Since the two outgoing 
nucleons are free, we have no particular spin selection rule in the de-excitation process for this case and thus all 
partial waves can contribute. As can be seen from Fig. |l^, the NP states become more important for the QF decay, 
and the overall energy shift more or less disappears. For the 2p breakup, the AA^ NN transition potential Vata 
provides not only a LO but also a TR component. The energy shift of the UP states is compensated by a relative 
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large enhancement of the NP states which results from TR excitation. Moreover, the p-wave rescattering interferes 
with the quasi-elastic scattering which has no intermediate AA^ configuration, hence the effects of Vaa are partially 
smeared out. 

To complete the discussion about the influence of the AA^ interaction, we show in Fig. |l^ the contributions of 
three characteristical AA^ partial waves to the exclusive ^H(p, 7i7r+)^H cross section. All selected partial waves have 
unnatural parity and therefore couple strongly to the pion. The spectrum is absolutely dominated by the "'iS'2 partial 
wave of the AiV system. Because the A and the nucleon have relative angular momentum zero in this case, there 
is no centrifugal barrier in the potential and the attraction becomes largest. That accounts for the strong energy 
shift in this partial wave. The partial waves with higher angular momentum L > 1 are subjected to a less attractive 
potential due to the centrifugal barrier and therefore do not exhibit such a lowering in the excitation energy. As can 
be seen from Fig. [l^, the ^P^ partial wave is only enhanced but not shifted due to Vaa- For the ^Dq partial wave, 
the potential has a minor effect and even gets repulsive. Compared to the higher partial waves are by far less 
important which results in the surviving of the downwards energy shift in the full spectrum. 

We conclude that neither the QF nor the BU process shows the same sensitivity to the LO channel as the CP 
process. Therefore only the exclusive ^H(p, n7r"'")^H spectrum clearly exhibits a lowering of the A excitation energy. 
This fact makes the coherent pion production most suitable to examine the effects of the AA^ interaction. 



D. Exclusive differential cross section for coherent pion production 



In this section, we present our results for the angular distribution of the coherently produced pions. In Fig. 14, the 
triple differential cross section a / dEndi^nd^-K of the exclusive ■^H(p, n7r+)'^H reaction is plotted as a function of 0^^ 
which is the angle between the outgoing tt"*" and the three-momentum transfer q. The energy transfer was chosen to be 
i^iab = 300 MeV and the neutron scattering angle is 9n = 0°. One recognizes that the angular distribution is strongly 
forward peaked. This means that most of the coherent pions are emitted into the direction of the three-momentum 
transfer. 

The contributions of the spin-longitudinal (LO) and the spin-transverse (TR) channel to the differential cross 
section are shown separately in Fig. |l^. One can understand some of the main features of these angular distributions 
by analyzing the spin structure of the transition operators. 

For the LO channel, the product of the de-excitation and excitation operators turns out to be 

(S* •_p;)(S"f • g) = • g • (pi X q). (43) 

The first term on the r.h.s. of Eq. (^) gives rise to a factor gp^ cos 9,^ in the transition amplitude and thus to a factor 
cos^ Ott in the cross section. This factor has a maximum at 0^ = 0° and thus partially explains the strongly forward 
peaked LO contribution in Fig. uB. There are, however, additional angular dependent factors such as the kinematical 



phase space factor and the overlap integral with the outgoing pion wave in Eq. (|39[). Those factors become larger 
as 0^ gets smaller and therefore pronounce the forward peaking of the angular distribution even more. All reaction 
events where the spin of the deuteron is flipped are produced by the second term on the r.h.s. of Eq. (|4^). At small 
scattering angles this cross section contribution is proportional to sin^ Oyr and thus vanishes at = 0°. 

We would like to mention that the angular distribution of the LO component is very similar to that of the pion 
elastic scattering |l|j2l[. In fact, one may view the coherent pion production as a kind of elastic scattering process, in 
which an initially off-mass-shell pion with the momentum q is converted into an on-mass-shell pion by the multiple 
scattering in the deuteron. This conversion process is possible because the deuteron as a whole can provide the extra 
recoil momentum needed to put the pion on its mass shell. 

For the TR channel, the de-excitation and excitation operators yield a spin structure 

_^ _ 2 i 

(S" •p;)(5^ X q) = -p; X q--{a X p^) x q. (44) 

As before, the first term on the r.h.s. of Eq. ( |4^ ) is connected with non-spin-flip events while the second term 
induces spin-flips of the deuteron. The angular distribution of the non-spin-flip part is proportional to |pi x g p = 
(^i^gsin^^)^. This factor vanishes for 0^ = 0° and peaks for 6*^ = 90°. The additional angular dependent factors 
as discussed before lead to a more forward peaked distribution. As can be seen from Fig. |l5|, the TR non-spin-flip 
component has its maximum at about 9-,^ ~ 40° and is zero at = 0°. The full TR angular distribution, however, 
has a different shape due to the spin-flip contributions to the cross section. This demonstrates the importance of 
the second term on the r.h.s. of Eq. (^4|) which does not vanish but rather reaches its maximum at = 0°. The 
overall result is quite a flat angular distribution at small pion angles and a fall off at higher angles which is less steep 
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as compared to the LO component. This behavior of the TR component is very similar to the observed angular 
distribution in pion photoproduction (7,7r) reactions on the deuteron, see e.g. Ref. [^ . 



VI. SUMMARY AND CONCLUSIONS 

In summary we have shown that (p, n) charge exchange reactions on a deuteron target provide an excellent tool 
to investigate the AA^ interaction. They probe both the LO and the TR response function in a kinematical region 
which is inaccessible to real pion or photon reactions. We have presented a coupled channel approach which allows 
the theoretical treatment of the interacting AN system. By assuming that the A resonance is excited dominantly in 
direct interaction with the projectile, we modified the coupled equations and showed how to solve for the correlated 
A7V wave function in a very efficient way with the Lanczos method. For the AiV potential, a meson exchange model 
was adopted which includes tt, p, w, and a exchange currents. 

Results of numerical analysis have been shown for inclusive and exclusive cross sections in the A energy region. In 
the coherent pion production ^II(p, n7r"'")^II, a strong shift of the A peak position is observed. We have shown that 
this lowering in excitation energy is due to the strongly attractive correlations in the LO spin-isospin channel. This 
attraction emerges mainly from the energy dependent tt exchange potential. Furthermore, we calculated the angular 
distribution of the coherent pion component and found it to be strongly forward (in the direction of the momentum 
transfer) peaked. 

The theoretical calculation for the inclusive ^H(p, n) cross section is in fairly good agreement with experimental 
data. The model describes both the A peak and the quasi-elastic peak very well. In the dip-region, theory and 
measurement agree only in the LO channel while experimental data are underestimated in the TR channel. The 
enhancement in the TR channel is most probably due to rescattering effects (i.e. purely mesonic exchange currents). 
If such effects could explain the shortcomings of the present model, they would also demonstrate the validity limits 
of the impulse approximation. Therefore rescattering effects are of special interest and need further investigation. 
In addition, our model may be improved by including higher nucleon (N*) resonances in the model space. Both 
rescattering effects and N* resonances will be subject of future studies. 
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APPENDIX A: THE EFFECTIVE PROJECTILE - TARGET-NUCLEON INTERACTION 

In the present paper, we approximate the spin-isospin part of t^N^NN in the NN cm. frame by 

tNN,NNis, t) = [a (cti • q){aj • g) + /? (ct^ x q){aj x q)]f^ ■ fj, (Al) 

where the unit vector q is connected to the initial and final nucleon momenta k and k' by q = it — k' . The coefficients 
a and P are functions of the Mandelstam variables s and t and describe the strength of LO and TR spin excitations, 
respectively. They are determined from experimental pn — > np scattering data |^ which can be fitted with 

where c„ = 208 MeV fm^, cp = 178 MeV fm^, A^i = 148 MeV, A„2 = 460 MeV, and Ap = 342 MeV. 
For the tNN,NA transition operator we assume in complete analogy the following form, 

tNN.NA{s,t) = [7 (a, • q){S} ■q)+-f'{a,X q){S] X q)]n ■ f} . (A3) 

Again, the coefficients 7 and 7' are adjusted to fit experimental data, here of the pN — > nA reactions 03], yielding 
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where = 488 MeV fm^ and A-y = 650 MeV. The fact that data are explained with 7 = 7' means that the ratio 
LO/TR is 1/2. 

The assumed operators are essentially of (5~function type; the only momentum dependence comes from the vertex 
form factors. In spite of their simplicity, tjs[N,NN and tMN,NA can reproduce not only the cross sections but also the 
spin observables from the reactions with the proton target pOl . 



APPENDIX B: EVALUATION OF MATRIX ELEMENTS 



All matrix elements which appear in this work are evaluated by using partial wave expansions and standard tensor 
operator techniques, see e.g. Ref. In this appendix, we discuss in some detail the calculation of the radial source 
function and of the matrix elements for the AA^ interaction. The matrix elements of Sec. IV are not given explicitly 
but can be calculated rather easily within the presented scheme. 



1. Explicit formula for the source functions 

For the AA'^ system, the radial source function is obtained from Eq. ( p^ ) by inversion 

PSLjMjir) = r {{SL)JM.j | pa) = r {{SL)JM.j | p | (Bl) 
where the deuteron wave function is given as 

1*'^)= E \K{r)\{lld)lMd). (B2) 

The hadronic transition operator p was defined in Eq. (^). If we apply the specific form of tj\iN,NA as given in Eq. 
(A3) and use the tensor operator notation, p may be rewritten as 

p = (IF) J2 PimAr) {SlY^'^ YUr) + [l ^ 2), (B3) 

where 

Pw(r)=47r7(.,t)(-r''+^(i ^^^Iw S \\\) ^' 3i{\qr)YC„M) (B4) 
and (IF) denotes the isospin factor 

(IF) - (A7V |(n I To • f t I p) I d) = 5iThMr^- (B5) 



Insertion of Eqs. (g|,g|) into Eq. (|BT]) yields 

PSLJMj{r) = (IF) ^ pimu{r) ^ (t>iA-r) 

Imv Ifi 

X {{SL)JMj I (Slt^Yi^if) I (l/d)lMd) + (1 ^ 2). (B6) 

The calculation of the angular momentum matrix element is straight forward, but requires some lengthy spin algebra. 
The final result is 

{{SL)JMj I (5j)Wll™(f) I {lld)lMd) = 

^ . JLild (-)^+l(<551 + ^5S2) E {-)^^'■'-^''\'^J^ + 1) 

JiMi 

I 1 ,h \( J, I J \ ( I Id L\l L ,h l\l S I l\ , . 

m Md -M, J \M, V ~Mj ) \ Q Q Q ) \ I I Id ] \ .h J L ] ' ^"^'^ 

For the radial source function of the NN system, an equation analogous to ( |B^ ) holds with 5^ replaced by a, an 
isospin factor (IF) = \/2, and a transition strength a{s,t) resp. P{s,t) instead of 7(5, t). 
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2. Matrix elements of the AA'^ potential 



In configuration space, the AN potentials Van^an (direct) and Van^na (exchange) as well as the transition 
potential Van~^nn exhibit the following spin-isospin structure, 



Here, a"" and r are the appropriate spin and isospin operators with a, a! £ {N, A} and 



(B8) 



(B9) 



is the usual spin tensor operator. In terms of the spin operators used in Sec. [II, we have ct^^ — t?, = S'^ and 
= thus the reduced matrix elements are given by 



-•NN 



\2 



-•NA 



b=2, (I II a^'^ II I) =2715. 



(BIO) 



Equivalent equations hold for the isospin operators r °° . 

The functions T of Eq. ( p38| ) are the Fourier transformations of the non-spin, spin-central and spin-tensor part of 
the interaction, i.e. 



27r2 
1 

2^ 

-1 



dk eUkr) [l^j^fc°,,,(fco, fc) + 2V;™ ,,(fco, k)] , 



•^^W = 2^ / ^'^2(fcr) [Vl:^,,Aka.k) - Vll,,,{k^, k)] 



(Blla) 
(Bllb) 
(Bllc) 



For the potentials given in Sec. [II, these expressions may easily be found analytically. Furthermore, we need the 
following spin matrix elements which are calculated with standard techniques. 



{[{sa'Sb-)S'L']JMj\ ar' -(ff I [{saSb)SL]JMj) = 

I \Sa+SU+S ! X i S 

{[{Sa'Si,')S'L']JMj\ S,%''''"''{f) I [{.SaSt)SL]JMj) - 



{Sa' II \\Sa){sb' Ik""' \\sb), 





L' 






S 




(o 





o) 




L' 


5} 



Sa' Sb' S' 
Sa Sb S 
1 1 2 



(B12) 



(B13) 



From Eqs. (BIO) - (B13) one can build up the complete expressions for the matrix elements Vnn'ir) as defined in Eq. 
m, i.e. 



Vnn'ir) - {[{Sa'Sb')S'L']JMj \ {{Ta'U')!! \ Vab^a'b' (ko , \ {TaU)!!) \ [{SaSb) S L]J Mj) . 



(B14) 
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TABLE I. Parameters used in the meson exchange model for the AiV interaction. 
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FIG. 1. Impulse approximation diagrams for the ^H(p, n) reaction. Here (a) shows the nucleon excitation and (b) the A 
excitation. 
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FIG. 2. Direct term (a) and exchange term (b) of the AN potential Vaa- The mesons taken into account are the pion (tt), 
the rho (p), the omega (w), and the sigma (a). 
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FIG. 3. Reaction mechanisms of the different contributions to the inclusive cross section in our analysis. Only the lowest-order 
diagrams arc shown. They represent: (a) quasi-elastic scattering, (b) p-wave rescattering, (c) coherent pion production, (d) 
quasifree A decay. 
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FIG. 4. Differential cross section for tlie reaction ^H(p, n) A"'"^ at Tp — 790 MeV and at scatteringangles 6n ~ 0°,7.5°, and 
15°, respectively. The theoretical calculation uses the parameterization of tNN,NA as given in App. M. Experimental data are 
from Ref. 
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FIG. 5. Zero degree neutron spectrum for the reaction ■^H(p,n) at Tp = 790 MeV. The theoretical calculation (solid line) 
includes contributions from quasifree A decay (QF, dashed), coherent pi on p roduction (CP, dotted), and two-nucleon breakup 
of the deuteron (BU, dashed-dotted). Experimental data are from Ref. |32| . 
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FIG. 6. Full calculation with inclusion of baryonic exchange currents (solid line) in comparison with the "spectator approx- 
imation" result where Vaa = VivA ~ (dashed line). Experimental data are from Ref. ]32|. 




FIG. 7. Separation of the inclusive cross section into spin-longitudinal (LO, ~ ct ■ g) and spin-transverse components (TR, 
' a X q). Experimental data are from Ref. [BSl. 
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FIG. 8. Neutron spectra for the reaction ^H(p, n) at Tp = 790 MeV at scattering angles of On 
respectively. Experimental data are from Rcf. |32[|. 
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FIG. 9. Effects on coherent pion production ^H(p, n7r+)^H resulting from the direct (dir) and the exchange (ex) contribution 
of the TT meson to Vaa (solid lines). The line labeled (dir+ex) shows the sum of direct and exchange contributions of the vr. 
Note the apparent shift as compared to the dashed line which represents the spectator approximation Vaa = 0. 
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FIG. 10. Effects on coherent pion production ^H(p, n7r"'")^H resulting from the 7r, p, a», and a meson contributions to Vaa 
(solid lines). The line labeled {it + p + u> + a) represents the full model calculation. The dashed line shows the spectator 
approximation Vaa = 0. 
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FIG. 11. Effects of Vaa on the exclusive cross section for coherent pion production (CP), quasifree A decay (QF), and 
deuteron breakup (BU). Calculations with (solid lines) and without Vaa (dashed lines). Note the different scalings of the 
vertical axis. 
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FIG. 12. Decomposition of the exclusive cross sections into contributions from unnatural parity states (UP) and natural 
parity states (NP). Calculations with (solid lines) and without Vaa (dashed lines). Note the different scalings of the vertical 
axis. 
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FIG. 13. Multipole decomposition of the zero-degree spectrum for the exclusive ^H(p, n7r"'")^H reaction. Calculations with 
(solid lines) and without Vaa (dashed lines) for three characteristical partial waves of the AN system are shown. Note the 
different scalings of the vertical axis. 
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FIG. 14. Triple differential cross section for the ^H(p, n7r+)^H reaction at wiab = 300 MeV and On = 0°. The cross section is 
shown as a function of On which is the angle between the outgoing tt"*" and the momentum transfer q on the deuteron target. 
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FIG. 15. Triple differential cross section for the ^H(p, n7r"'")'^H reaction in the spin-longitudinal (LO) and the spin-transversal 
(TR) channel. Mi (Mf) refers to the spin-projection of the deuteron in the initial (final) state, hence the dashed curves labeled 
Mi = Mf show a calculation where spin flips of the deuteron have been excluded. 
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